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O { - We calculate the QED and QCD radiative corrections to the charged lepton 

energy distributions in the dominant semileptonic decays of the top quark t — > 
6VF + — ► b{i + V() (I = e, fi, t) in the standard model (SM), and for the decay 
t — > bH + — > b(r + v T ) in an extension of the SM having a charged Higgs boson 

o 



next-to-leading logarithmic approximations, but the QED corrections are considered 
in the leading logarithmic approximation only. These corrections are numerically 



important for precisely testing the universality of the charged current weak inter- 



actions in i-quark decays. As the r + leptons arising from the decays W + — ► r H 



and H + — > t + u t are predominantly left- and right-polarised, respectively, influenc- 

5-H ' 



ing the energy distributions of the decay products in the subsequent decays of the 
r + , we work out the effect of the radiative corrections on such distributions in the 
dominant (one-charged prong) decay channels r + — > tt + u t , p + v T , a\v T and l + U£i> T . 
The inclusive tt + energy spectra in the decay chains t — ► b(W + ,H + ) — ► 6(r + i/ T ) — ► 
b(n + v T v T + X) are calculated, which can help in searching for the induced effects 
at the Tevatron and the LHC. 
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I. INTRODUCTION 



Top quark is now firmly established by the experiments CDF and DO at the pp collider 
Tevatron at Fermilab, with m t = 173.1 ±1.4 GeV, decaying dominantly through the mode 
t — > bW + — > b(£ + is£, qq') [1]. At the Large Hadron Collider (LHC), expected to be operational 
shortly, one expects a cross section a(pp — > ttX) ~ l(nb) for the LHC centre of mass energy 
of 14 TeV [2]. With the nominal LHC luminosity of 10 33 (cm)~ 2 (sec) _1 , one expects a ti 
pair produced per second. The ti production cross section for the 10 TeV run of the LHC is 
estimated as about 0.4 nb 2], still large enough to undertake dedicated top quark physics. 
Thus, LHC is potentially a top factory, which will allow to carry out precision tests of the SM 
and enhance the sensitivity of beyond-the-SM effects in the top quark sector. Anticipating 
this, a lot of theoretical work has gone into firming up the cross sections for the tt-pair and the 
single-top production at the Tevatron and the LHC, undertaken in the form of higher order 
QCD corrections 3], 0,0, (|. Improved theoretical calculations of the top quark decay width 
and distributions started a long time ago. The leading order perturbative QCD corrections 
to the lepton energy spectrum in the decays t — > bW + — > b(£ + U£) were calculated some thirty 



0. 



years ago |Y|J. Subsequent theoretica 
0(a s ) corrections were published in 



work leading to anal ytic derivations implementing the 
9| and corrected in [lOj . The order a s contribution to 



the top quark decay width dominates the radiative corrections (typically -8.5%). The 0(a) 



12 1, the finite W-width effect 



electroweak corrections contribute typically +1.55% [ll|, 

1.56%) almost cancels the electroweak correction [3]. The next-to-leading order (N 
QCD corrections in a s (i.e., a 2 s ) were computed as an expansion in (Mw/m t ) 2 in 
These results were confirmed later by an independent analytic calculation in 16j, [17[ , and 
contribute about -2.25% to the top quark decay width. 

Our main concern in this paper are the lepton energy distributions from the decays 
t — > bW + — > bijt^vi) (for £ + = e + ,/i + ,r + ), which are modified from their respective Born- 
level distributions in a way specific for each charged lepton due to the QED corrections. 
These (QED and QCD) radiative effects have to be taken into account to test the universality 
of charged current weak interactions in the top quark sector. Another process which breaks 
the charged lepton universality in the decays t — > bl + v t is induced by charged Higgses 
(for m,H+ < m t — m^) in the intermediate state, t — > bH + — > 6(£ + ^), which is expected 
to influence mainly the final state br + v T due to the H + t + vi couplings. The leading order 
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in a s corrections to the polarized top quark decay into H + b have been calculated in 18]. 
We study the effects of the radiative corrections on the r + -energy distribution in the decay 
t -> bH + -> br+v T . 

Radiative (QED and QCD) corrections in the top quark decays, such as t — > hW + — > 
b£ + v#, with £ + = e + ,fi + ,r + , involve large logarithms due to the large fermion mass ratios. 
For example, in the leading logarithmic approximation (LLA), one encounters the logarith- 
mic terms 



2 x / m 2 



L e = lngj«25.4, L M = ln^J«14.8, (1) 
L T = In ( 4) « 9.1, L b = In ( » 7.4, 

in the partial decay widths. Hence, in the LLA, radiative corrections to the partial widths 
lead to typically large effects 

ry ry ry ry 

-L e ^6.2%, -L^3.6%, -L T &2.1%, —L b ^23%. (2) 

7T 7T 7T 7T 

They are included together with the non-logarithmic terms in the estimates undertaken in the 
fixed order (in a or a s ) calculations. However, to get perturbatively reliable results, all terms 

/ 2 \ n 

of the type (^) n In f ^ J in the decay t — > 6e + z/ e , for example, have to be summed up (the re- 

summed leading log approximation LLA), as well as (^) n In f ^| j (the next-to-leading log 
approximation NLLA). Using the well-studied case of the QED radiative corrections to the 
purely leptonic decays — > z/ M e - z/ e , we show that the structure function (SF) approach 19, 



201 ] (based on the factorisation hypotheses [2l|]) is the appropriate framework to resum such 



terms, enabling us to derive the electron energy spectrum with the radiative corrections taken 
into account to all orders of the large logarithms. As a warm-up exercise, and also to set our 
notations, we reproduce the well-known results for the QED corrections to the muon decay 



22 



23 



24 



251 ] and generalise it to all orders of perturbation theory by summing 



up the leading logs (^liifm^/m^))" (see Section EJ) . In this context, we also discuss the 
polarised muon decay case. The SF approach is applied next to the semileptonic decays of the 
top quark t — > bW + — > b{[ + v^), where the QCD and QED radiative corrections to the Dalitz 
(double differential) and inclusive lepton energy distributions are worked out. In this, the 
QCD-corrected energy distributions are derived in the re-summed leading logarithmic and 
next-to-leading logarithmic approximations, but the QED corrections to these distributions 
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are calculated in the leading logarithmic approximation only. This is discussed in detail in 
Section IVl 

In many extension of the SM, the Higgs sector of the SM is enlarged, typically by adding 
an extra doublet of complex Higgs fields. After spontaneous symmetry breaking, the two 
scalar Higgs doublets $1 and $2 yield three physical neutral Higgs bosons (h, H, A) and a 
pair of charged Higgs bosons (H ). If m#+ < mt — ra^, one expects measurable effects in 
the top quark decay width and decay distributions due to the i7 ± -propagator contributions, 
which are potentially large in the decay chain t — > bH + — > b(r + u T ). The two parameters 
which determine the branching ratio for this decay are and the quantity called tan /3, 
defined as tan/5 = V2/V1, where v% and t>2 is the vacuum expectation value of $1 and $2, 
respectively. Of particular interest is the parameter space with large tan (3 (say, tan (3 > 20) 
and m,H± < 150 GeV. This mass range is already excluded (for almost the entire tan/3 
values of interest) in the so-called two-Higgs-doublet-models 2HDM due to the lower bound 
on m#± of 295 (230) GeV at the 95%(99%) CL. from the experimental measurements of 
the branching ratio B(B — ► X s ^) [26], and the order a 2 s estimates of this quantity in the 
SM [271 ] . However, this bound applies only to the 2HDM of type II, in which the Higgs 
doublets $1 and $ 2 couple only to the right-handed down-type fermions (dm, £m) an d the 
up- type fermions (um, Um) , respectively. In the minimal supersymmetric standard model 
(MSSM), one has a type II 2HDM sector in addition to the supersymmetric particles, in 
particular the charginos, stops and gluinos. Their contributions could, in principle, cancel 
that of the charged Higgs bosons in the B — > X s + 7 decay rate. Hence, the 2HDM-specific 
constraint on m#± from B(B — > X s ^) is not applicable in the MSSM. In our opinion, the 
natural embedding of the extra Higgs doublet is in a supersymmetric theory, and hence 
we will ignore the lower bound on m#± from B(B — > X s j). A model- independent lower 
bound on m H ± exists from the non-observation of the charged Higgs pair production at 
LEPII, yielding m#± > 79.3 GeV at 95% C.L. 26], which we shall use in our numerical 
analysis. Thus, a charged Higgs having a mass in the range 80 GeV < vnu± < 160 GeV is 
a logical possibility and its effects should be searched for in the decays t — > bH + — > r + v T . 



A beginning along these lines has already been madejit the Tevatron 
definitive search will be carried out only at the LHC 



31 



28 



29J, laol 



29, 30], but a 



321 ] . We work out the effects of 



the radiative corrections to the lepton energy spectra in the decays t — ► bH + — ► b(r + v T ) in 
Section ED 
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The r + leptons arising from the decays W + — > t + v t and H + — > r + z/ r are predomi- 
nantly left- and right-polarised, respectively. Polarisation of the influences the energy 
distributions in the subsequent decays of the t ± . Strategies to enhance the if ^induced 



effects in the decay t — > bW + — > 6(r + z/ r ), based on t 



33 



re polarisation of the r + have been 



34 



35 



, 136 



371 ] . We work out the ef- 



discussed at length in the existing literature 
feet of the radiative corrections on such distributions in the dominant (one-charged prong) 
decay channels r + — > 7r + z/ r , p + z/ r , a\v T and £ + z/^z/ T . To implement this, we again use the 
SF approach 38|. In particular, the inclusive ir + energy spectrum in the decay chain 
t — > b(W + ,H + ) — > 6(r + z/ r ) — > b(7T + u T u T + X), and likewise for the decay chain of the i 
quark, can be used to search for the induced effects of the H ± at the LHC and Tevatron. 
Details are given in Section I VI I and in Appendix A. 

To get the relative normalisation of the decay width t — * bH + with respect to the SM 
decay width t bW + , one has to take into account the loop corrections (quantum soft 
SUSY-breaking effects). These quantum effects on t — > bH + have been worked out in the 
context of the minimal supersymmetric standard model MSSM in a number of detailed 
studies (see, for example 39|, |40|]), and the bulk of them can be implemented by modifying 



the 6-quark mass, 777,™ rrected = mb j (\ _|_ /\ b y The specific values of depend on the super- 



symmetric mass spectrum, and can be calculated using FeynHiggs 41], given this spectrum. 
The influence of these corrections on the branching ratio for the decay t — > bH + have been 
recently updated in [42], predicting BR{t — ► bH + ) > 0.1 for < 110 GeV in the large- 
tan p region (tan (3 > 40). We shall pick a point in the (tan (3 — m^+) plane from this study, 
allowed by all current searches, for the sake of illustration. We summarise our results in 
Section EH 



II. MUON DECAY: A WARM-UP EXERCISE 



We start by discussing the electron energy spectrum in \x — > eVgV^ decay. In the Born 
approximation, this spectrum is given by the following formula 431 ]: 



B 



dY 

dx 



6r, 



2x 2 (I - x) 



9 



px 2 (3 — Ax) 



(3) 
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where x = 2E e /m^ is the energy fraction of final electron, p is the well-known Michel 
parameter 



44] and T is the total decay width: 

r„ 



1927T 3 ' 



(4) 



where Gp is the Fermi coupling constant. Using the SF approach 19|, |20), we can derive the 
electron-energy spectrum with the radiative corrections taken into account to all orders of 
the large logarithm: 

l 

(5) 



dT 



RC 



dx 



dy ^ fx A dT B / H a T ^ . A 



a 



Z7T 



L = In ( — I 

mi 



10, 



where ^ is the electron spectrum in the Born approximation ([3]) which is considered as the 
hard sub-process. D(x,/3) is the so-called structure function, which describes the virtual 
and real photon emission in the leading logarithmic approximation and has the form 381 ] : 



D (x, (5) = 6 (1 - x) + (3P {1) (x) + -ft 2 P {2) {x) + • • • . 



(6) 



The quantities (x) are the kernels of the evolution equations which are defined by the 
following relations: 

1 



P (1) (x) 

p (n) fx) 



1 — x 



lim 



i±^!^(l-x-A)+ (2 In (A) + 
1 — x V 2 



<Hl-x) 



, (7) 



dy 

y 



p(l) > y ) p(n-l) 



The structure function D (x, (3) defined in this way automatically satisfies the Kinoshita- 
Lee-Nauenberg (KLN) theorem 4a, |46j on the cancellation of the mass singularities in the 
total decay width 



dxD (x,P) = 1. 



(8) 



There also exists a smoothed form for the structure function D (x,/3): 
D (z, (3) = 2/3(1- zf- 1 (l + lpj-p(l + z ) + {(3 2 ) 



(9) 
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which sums radiative corrections in all orders of perturbation theory which are enhanced by 
the large logarithmic factor L (in f3) and is more convenient for numerical evaluation. 

The quantity K (x) in (jHJ) is the so-called i^-factor which takes into account the con- 
tributions of the radiative corrections which are not enhanced by the large logarithms and 
have rather complicated form (see 22J or 47J], §147). We note that, contrary to the singular 
behaviour (~ In (1 — x)) of K (x) in the limit as x — > 1, the quantity 



(10) 



has a finite limit as x — > 1 481 ]. 



Thus, applying the general form of the corrected spectrum ([5j), we obtain the following 
form of the electron energy spectrum in the leading logarithmic approximation (LLA): 



1 dT 

6T dx 



2x 2 



1 — x — -p (3 — Ax) 



+ 



aL 
2^ 



4Fi (x) 
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pF 2 (x) 



(11) 



where the functions F\ } 2 (x) are the results of the application of the structure function to 
the spectrum in the Born approximation: 



F 1 (x) 



F 2 (x) 



2x 2 (l-x)ln 
l 

/% 2 (3-4 y )P« 

J y 



(- — -J + - (1 - x) (1 + ix - 8x 2 ) 
\ x J 6 ' 



(12) 



= 2x 2 (3-4a;)ln 
which satisfy the following property: 



x 



x 



16 o 2 1 

3 o 



(13) 



dxFi ^ (x) = 0. 



(14) 



This is a specific form of the general KLN theorem 45|, 146 ]. 

In concluding this section, we give the double differential distribution for the case of the 
polarised muon decay with the radiative corrections in LLA (here we put p = 3/4): 



dT 



Tdxd cos 9 



= x 2 [3-2x- P At (l - 2x) cos6] 
+ ^[F 3 (x)-P fl coseF 4 (x)} 1 



(15) 
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with and 9 being the degree of muon polarisation and the angle between the muon polar- 
isation vector and the electron momentum (in the rest frame of the muon). The functions 

F 3 = 2 (6F 1 - F 2 ) , F A = 2 (-2F 1 + F 2 ) , (16) 

have the explicit expressions: 

F 3 (x) = Ax 2 (3 - 2x) In ^— - + | + Ax - 8x 2 + ^-x 3 , 

1 — x 1 

F 4 (x) = 4x 2 (1 - 2a;) In Ax 2 + 8a; 3 . 

x 6 

III. TOP QUARK DECAYS t -» 6(W + ,#+) IN THE BORN APPROXIMATION 

Top-quark decays within the Standard Model are completely dominated by the mode 

t^b + W + , (17) 

due to Vtb = 1 to a very high accuracy. In beyond-the-SM theories with an extended Higgs 
sector, if allowed kinematically, one may also have the decay mode 

t^b + H + (18) 

where H + is the charged Higgs boson, which we will consider within the MSSM. The relevant 



part of the interaction Lagrangian is 49]: 



Cl = 2wk^ VtbH+ [Ut {pt) {A (1 + 7s) + B (1 ~ 7s)} Ub {Pb)] + 

+ ^f^ H+ K M (1 - 7s) u, (p,)] , (19) 

where A, B and C are model-dependent parameters which depend on the fermion masses 
and tan j3: 

A = m t cot/3, B = nit, tan j3, C = m T tan (3. (20) 



The decay widths of processes ( ITT)) and (ITS)) in the Born approximation are well known 49] : 



<T \i / m b M\ 



Tt-*vw — -7T. — 7J2 ( 1) — |) — T 

647rM^mi \ 

g 2 .1/ m 2 



(m 2 + m 2 ) + (m 2 - m 2 ) 2 - 2M, 



, (21) 



Ft-^bH — — 772 A 2 1, 2" ' X 



x 



[(m 2 cot 2 /3 + m 2 tan 2 (3) (m 2 + m\ - M 2 H ) - Am 2 m 2 b ] , (22) 
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where A (x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz is the triangle function. The total top 
quark decay width then reads as: 

r t °* = ^t-^bW + ^t-^bH- (23) 

We now discuss the total top quark decay width including the radiative corrections. In the 
total decay width the contribution of the QED corrections containing the large logarithms L 



is cancelled (see ((8 
were calculated in 



50 



he non-enhanced QED corrections are small. The QCD corrections 



5]J and have the form: 



-ptot j^Born+QCD , -pBorn+QCD ( r )A\ 

1 t,RC — 1 t~>bW ~r 1 t^bH > 

v Born+QCD _ p /, , f x , _ / \ 4lT 2 \ 
1 t->b(W,H) ~ 1 t-*b(W,H) I 1 + JW,H) , JW,H — — I — j ■ 

IV. THE TOP QUARK DECAY t -> bW + -> b(£+v e ) IN THE BORN 

APPROXIMATION 

The formalism illustrated in Section [IT] can be used to discuss the inclusive semileptonic 
decays of the charm, beauty and top quarks. However, the decay distributions from the 
the charm and beauty hadrons have in addition important non-perturbative effects, which 
usually are modelled in terms of the shape functions. In the case of the top quark decay, since 
the top quark lifetime is much shorter than the typical strong interaction time, the decay 
dynamics is controlled by perturbation theory. Thus, incorporating the (QED and QCD) 
perturbative corrections, one has precise theoretical predictions for the energy spectra of 
the decay products to be confronted with data. We start by working out the charged lepton 
energy spectra in the decays t — > bW + — > 6(£ + z/^), where £ + = e + , r + . To that end, let 
us consider the dominant decay in the SM (see Fig. [Tj, a.)): 

t (ft) -> b (p 6 ) + W + (q) -> b ( Pb ) + (£+ (p e ) + v, (p„)) , (25) 

and to be specific, we concentrate on the case with l^v^ = e + v e . The matrix element of this 
process in the Born approximation is given by: 

t ->bW+->b(e+v e ) _ , 9 2 Vtb 1 f q^qu\ 



4V^9 2 -^r M 2 W/ 
x [u b ( Pb ) 7 M (1 + 75) Ut (pt)} [u e (Pe) -f (1 - 7s) u Ve { Pv ]\ , (26) 
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a.) 




b.) 



Fig. 1: Lowest order Feynman diagrams describing the semileptonic decays of the top quark a) 
SM, mediated by W + , b) BSM, mediated by H+. 



where g 2 



47TQ 

sin 2 6w 



SnM^G p I 'v2 is the electroweak coupling constant, 6w is the weak 



mixing angle, Myy is the W /± -boson mass, and V t b is an element of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix 52J, |53(. We note that the contribution of the second 
term in the parenthesis is proportional to the electron mass due to the conservation of the 
lepton current and can be omitted. The matrix element squared then reads as: 



M 



t^bW+^b(e+u e ) 



Born 



9 2 V t 



ib 



4V2(g 2 -M^; 



2 8 (PbPu) (ptPe) ■ 



(27) 



Let us introduce the following notation for the kinematic variables: 



Xfo ra t ' 
M w ' 



7 



2E e 



V 



2E V 



where Eb, E e and E v are the energies of the 6-quark, positron and neutrino in the t-quark 
rest frame, respectively. T w is the total decay width of the VT-boson. In terms of these 
variables, the various scalar products can be expressed as: 



2 (PbPe 

2 {ptPe 



mi [I 



0C Jj 



v) 



2 {pbPv 
2 (ptPu 



mi (I 



mrx 



7]), 2{p e p v ) = m 2 t (I + 7] - x b ) 
2 (ptPb) 



m\xb- 



Since the main contribution to this decay comes from the kinematic region where M^-boson 
is near its mass-shell we have to take into account its decay width. We use the Breit-Wigner 
form of the propagator: 
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Thus, the matrix element squared (127)) then reads 

I 1 T 7)1 

e 2 . (29) 



t_*6 W H-_>6(e+i/ e ) 
iW Born 



2 2(^ 2 ^) 2 x e (l-a; e -r / ) 



The phase space volume element with three-particle final state has the standard form: 

rf$ 3 = (27T)- 5 5 (ft - Pb - Pe - Pv ) = ^3^e^- (30) 

The kinematic restrictions are: 

< x e < 1 — rj, 

1 — X T < Xb < 1, 

1 — x e — t] < x u < 1 * 



1 - x ' 



and the 6-quark mass-shell condition fixes the cosine of the angle between the positron and 
the neutrino momenta directions C eu = cos (9 ev ) = £ e *Jf | , 

2 

C ej , = H (1 - x e - x u - rj) . (31) 

X^Xii 

On using the standard formulae for the decay width 

dT = — - — |M| 2 d$ 3 , (32) 

we obtain for the case of the unpolarised top quark decay t — > bW + — > b{i + vi) the decay 
width: 

dT t ^ + ~* Kl+Vl) _ xi (1 - xi - rj) = xx {xf ax - x t ) 



dx h d Xl (1 - ^ (1 + ?7 - Xf,)) 2 + 7 2 ( +7 2 



/2 (I" 

where y = 1 + rj — x b , x™ ax = 1 — rj and y Q = T t is the dimensional factor: 

r, = «I. (34) 

Now, we calculate the branching ratios of the decays considered above. The branching 
ratio of the decay t — > bW + — > 6(£ + z/^) is obtained from (133)) by dividing it by the total 
width of top quark Y\ ot (see Q23J): 

dx b d Xl Vl_i) 2 +7 2' * 
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Let us consider the electron energy spectrum. In the Born approximations it has the follow- 
ing expression: 

i 



where 



dBr® dBr t r bW+ ^ b{e+Ue) 



Born 



B t ■ x e (x r e 



dx, 



dBr 



t->bW+->b(e+9 e ) 
B 



doc doc g 



1— X e 
X e ) $ W (%e 



dy 



l-JL 

yo 



1 

7^ 



arctan 



) + 7 2 



arctan 



7 



C(r) + x) -1 

7 



(36) 



(37) 



A. QCD radiative corrections 



The inclusive electron energy spectrum including the lowest order QCD corrections is 



dBr 



t^bW+^b(e+v e ) 
Born+QCD 

doc p 



i / dr(°) + dv^ CD 



~ptot 
1 t,RC 



doc p 



(38) 



where r*°jj C is the radiatively corrected total decay width of top quark from (1241) . This 
expression is free from the 6-quark mass singularities, hence we can put r\ — 0, which yields: 



dr {1) 

al QCD 



2a, 
3vr 



dy 



+7 : 



: F W {x e ,y) 



(39) 



where the function F\y (x, y) is finite in the limit m& — > and has the form 101 ]: 



F w (x, 2/) = 2a; (1 — rr) 



C2+Li2(l)+Li2( |) + iw(i^ 



+ 



+ x 



y 

( 2 + Li 2 (?/) - Li 2 (x) - Li 2 ( - 



+ 



+ -In (1-2/) [-(3 + 2a;)+2 2 /(l + x) + 2 / 2 ] + 

+ - In C 1 - [a: (9 - 4s) - 2y (1 + a?) - y 2 ] + 
2 V x/ 

+ ^^ln(l-x) + ^(l-x) + 



(40) 



2 v ' 2- ■ ■ \x 

This formula is valid for x e < 1. For x e ~ 1, close to the boundary of the phase space, there 
are Sudakov Logarithms due to the limited phase space, and this result becomes unstable 
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Fig. 2: Kinematics depicting the application of the structure function method, which involves 
factorisation of the amplitude in the "hard sub-process" (filled circle) and the "long-distance" 
contributions (empty circle) taken into account by the convolution with the structure function 
DfaP) (see (02}. 



but remains integrable. The electron energy spectrum with the QCD corrections is given 
by: 

_t-»6W+-»&(e+f e ) 



dBr 



Born+QCD 



Ft 



dx P 



dy 



2a 

X e {xT X ~ X e ) ~ (x e , y) 

07T 



(41) 



B. QED radiative corrections in the leading logarithmic approximation 



To calculate the QED radiative corrections in the leading logarithmic approximation we 
will use the SF method which was illustrated in Section [III (see kinematic scheme in Fig. [2l 
a). 

The QED radiative corrected spectrum is: 



dBr 



t^bW+^b(e+9 e ) 
Born+QED 



V tot I ii 

1 t,RC j y< 



dv fx \ dT ^W+-,b(e+u e ) 



dy e 



(42) 



where the structure function D (x,/3 e ) was denned in (jBj). 

The first order QED radiative correction reads as (using the electron energy spectrum in 
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the Born approximation (|36j) ): 

jo t-+bW+-+b(e+i>e) 1 /• J / \ 7T^t^&W+^6(e+P e ) 

dBr QED LLA = X \ 1 f ^£p(l) J ^£ ] dT B 

dx e 2tt e r*^ c J y e \y e ) dy e 



2tt V 



i 

1 t,RC ^ i/e \ ye / 



t,RC 



where y™ ax = 1 — rj, and 
l 

IV) = J jP {1) {^)y{y max -y)^w{y) 



{x) <x(l — x) 



, . 1 — x\ 3 
+ 2 



+ x In (a;) + (1 — a;) 2 — - (l — x 2 ) 



i 

+ / d y^—^J^i^- l*w (V) - *w (x)] , (44) 



X 



dBr^ 

where $iy (x) is given in (|37|) . The contribution of the QCD correction — d ® CD from (|38l) is 
shown in Fig. [31 and is the same for I = e, fi, r. The contributions of the QED corrections 
is specific to the charged lepton e, /i, r and shown in Fig. [31 The input parameters used 
in this figure and subsequently are given in a table in the Appendix. In Fig. HI we show 
the electron energy spectrum in the Born approximation and compare it with the (QED + 
QCD) radiatively corrected ones. 

It is obvious from the foregoing that the QED radiative corrections break the lepton 
universality, encoded at the Lagrangian level for the decays t — > bW + —>■ b£ + i>£. It is also 
clear that the radiative corrections are not overall multiplicative renormalizations and they 
distort the Born level distributions in a non-trivial way. To quantify this, we plot the ratios 
R eT (x) and i?^ T (x), defined below, in Fig. [51 



R, 



er(x) 



Born+QCD+QED LLA 



\ L \ x — e / / Born+QCD+QEDLLA 
( V t^bW+^b{r+v T )f x _ x \\ 
p _ T / Born+QCD+QEDLLA 



Born+QCD+QEDLLA 
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Fig. 3: QCD and QED corrections to the lepton energy spectrum in the decays t — > bW + — > 

dBr 

b(e + u e ,T + i' T ). The solid curve is the QCD-correction term — J^ CD (i.e. second term from from 

(|4ip). the dashed and dotted curves are the QED-corrections from (|43p for 

the r + and e + in the final state, respectively 

As can be seen, the effect of the radiative corrections is very marked for the low- a; values 
of the lepton-energy spectra {x < 0.3) and it is non-neglible also near the end-point of the 
spectra (x > 0.7). This is numerically an important effect and in the precision tests of 
the SM in the top-quark sector, which we anticipate will be carried out at the LHC, it is 
mandatory to take the radiative distortions of the spectra into account. 

V. THE TOP QUARK DECAY t -> bH+ -» b(l+&i) IN THE BORN 

APPROXIMATION 

Let us consider now the top quark decay induced by a charged Higgs boson: 

t (Pt) - b (p b ) + H + (q) ^ b (p b ) + (£ + (p e ) + v, (p,)) (46) 
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Fig. 4: Lepton energy spectra from the decays t — > bW + — » b((. + ve) versus the lepton energy fraction 
x = x e ^ )T . The spectrum in the Born approximation (solid curve) is the same for £ + = e + ,/x + ,r + 
(see (f33|) ). The dotted curve is the e + -energy spectrum (dBr/dx e ) including the (QCD+QED) 
radiative corrections for the decay t — > bW + — ► b(ev e ) (i.e. the contributions from (I4ip plus the 
QED correction term from (|43p ). The dashed curve is the r + -energy spectrum (dBr /dx T ) including 
the (QCD+QED) radiative corrections for the decay t — > bW + — * b{Tv T ). 

where we will concentrate on £ + = r + (see Fig. [TJ b). Using the couplings from the La- 
grangian ffl~9]) we can write the matrix element of the process (|46p in the following form 

x [u t (p t ) {A (1-75)+ B(l + l5 )}u b (p b )}. (47) 



The model parameters A, B, C are given in (120]) . Squaring this matrix element yields 



M t^bH+^b{r+u T ) 



= (s »-j^)'+j^n, (BPl) (p *> a* • (48) 



Introducing the kinematic variables: 

i T = x„ = u = -4 = 1 4 n - x 6 , 

_Th_ 
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Fig. 5: Ratios of the lepton-energy spectra in the decays t — ► — > 6(£ + i^): R eT {x) = 

ir -'":; 1 "-' / J, "'t"°" (- ™). ™<j w = J *^-V' fl, '""^ < " ) ( d - hed 

curve), quantifying the leading order (QCD and QED) corrections to the lepton universality in 
semileptonic top quark decays. 

where E T and -E^ T are the energies of the final r lepton and the neutrino in the t-quark 
rest frame, respectively. Th is the total decay width of the charged Higgs boson, and the 
Breit-Wigner form of the propagator reads as 

1 111 



\ q 2 _ M 2 h + iMh T h \ z (g 2 - Ml) 1 + m*t% M%^_±y + ^ 
Thus the matrix element squared (148|) takes the form: 



(49) 



M t^bH+^b(r+u T ) 



= [9 V tb) T7Z T2 7-2- (50) 

M w +7 ^%o 

The decay width of the unpolarised top quark decay t — > bH + — > b(r + u T ) then takes the 
form: 

dv t^H^r^ T ) ^ {xmax _ ^ 



rf / v N2 ° J , (si) 



r * = oTTZi 5T7i (gVtb) T7A- I 52 ) 



2 u tt 3 V M§v J vy 7 M£' 
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The branching ratio of this decay is: 



dBr 



t^bH+^b(r+u T ) 



B 



doc jj doc *r 



ii 



1 - i 

3/0 



+ 7h 



1 1 



(53) 



where is the branching of the decay H + — > r+iv 49(: 



5f 

^H-*TV T 



32ttM^ M ' 



(54) 



3t? 2 M g 
32ttM2, 



(il4cot 2 /3 + Mjtan 2 /3) . 



For the numerical values of tan /3 that we entertain in this paper, the branching ratio = 1, 
to a very high accuracy. The dependence of the branching ratio of the decay t — ► M/ + on 
tan p is plotted in Fig. El We emphasize that in plotting this figure, radiative corrections 
coming from the supersymmetric sector are not included. They have been calculated in great 
detail in the literature, in particular for the MSSM scenario in 40], and can be effectively 
incorporated by replacing the 6-quark mass m;, in the Lagrangian for the decay t — > bH + by 
the SUSY-corrected mass m ,c° rrected = + AJ. The correction A& is a function of the 

supersymmetric parameters and, for given MSSM scenarios, this can be calculated using the 



FeynHiggs programme 



40| | . In particular, for large 



4=11 ] . which makes use of the results in 
values of tan/3 (say, tan (3 > 20)), the MSSM corrections increase the branching ratio for 
t — > 6if + significantly though this is numerically not important for tan/? = 22, which we 
use to numerically calculate the branching ratio for t — > bH + . We emphasize that in the 
analysis of data in the MSSM context, the branching ratio shown here in Fig. [6] has to be 
corrected to include the SUSY corrections. This, for example, can be seen in a particular 
MSSM scenario in a recent update [42], based on the version FeynHiggs v2.6.2. 

The lepton energy spectrum in the Born approximations has the following expression 

J x b B , , = B*B? ■ $ H (x T ) , (55) 



dBr 



t^bH+->b(T+u T ) 



13 



doc cIjQlj j- 



1— X T 
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Fig. 6: Lowest order branching ratio for the decay t — ► bH + as a function of tan/3 for M H + = 120 
GeV. 



where 



dy- 



y(i -y) 



y 2 



l — 2/o 



) +7 2 



7h 

+ (2y -l) In 



, 1 \ f x — Do 
arctan ( — + arctan 

IhJ V Vol 



JJo 



+ 



yo - x 



(56) 



A. QCD radiative corrections 



The leading order QCD corrections to the decay t — > bH + — > b{r + v T ) is calculated in 
a similar way as for the case of t — > bW + — ► b(r + u T ). The derivations for the Dalitz 
distribution dSr / dxf,dx T and the r-energy spectrum cLBr / cte T are given in Appendix [B] 
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B. QED radiative corrections in the leading logarithmic approximation 

To calculate the QED radiative corrections in the leading logarithmic approximation we 
will again use the structure function method, which gives: 

Born+QED _ & T \ U>V p ( fr ^ \ "i B ^7) 



dx T r™£ c J y \y J dy 

where the large QED logarithm now is (3 T (see (J3])). The first order QED radiative correction 
reads as 

j n t-*bH+-*b{l+v{) R 1 / \ r r t-^bH+-^b(l+ui) 

dxi 2vr ' r*°^ c J y \y J 



dy 



-(^-1)^1 



B "Jj pm (?) * hW 



t,RC 



( L i ~ 1) ndH^ (*0 - (58) 



where (x) is the Born spectrum defined in (I56I) and 



1 fl — X 

$ H (x){x+- + \n(x) + 2\n 



2 v 7 V x 



1 

.2 1 ^2 



+ / dy { y +x) ' [$ g (y) - $ g (s)] . (59) 

The contribution of the QED corrections ( 1581) is shown in Fig. [7] and compared with the 
QCD corrections from (1B23I) . In Fig. [8] we show the lepton energy spectrum in the Born 
approximation and compare it with the radiatively corrected one. Contrasting the r-energy 
spectra in this figure with the corresponding spectra in Fig. H] shows that the r-leptons from 
the decay t — ► bH + — > br + v T are distinctly more energetic. This feature is well known in 
the literature. We have calculated here the (QCD + QED) corrections to these spectra and 
checked their perturbative stability. 
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Fig. 7: QCD and QED corrections to the lepton energy spectrum in the decays t — > bH + — * bT + v T . 
The solid curve is the QCD corrections, i.e. the second term on the r.h.s. of the first line of 
Eq. (|B24jl divided by the total decay width T\ ot RC , and the dashed curve is the QED corrections 

dBr t^bH+ ^br+ V T 

QED ( to A from Eq. ([58]) for the r + in the final state. 

VI. TOP DECAY CHANNELS INVOLVING THE r LEPTON 

The radiatively corrected charged lepton energy spectra from the decays t — > bW + — > 
bi + vi (for £ + = e + ,/i + ,r + ) and t — > bH + — > b{r + v T ) presented here will be helpful in 
undertaking precision tests of the SM and in the searches for the if ^induced effects in the 
semileptonic decays of the top quark. Integrating these spectra from some experimental 
threshold lepton energy, the anticipated enhancement in the branching ratio for the t — > 
b{r + u T ) mode over the other two semileptonic modes t — > + u^, e + z/ e ) provides the 
experimental handle on the H ± searches. This is the strategy which is being used at the 
Tevatron, where searches have also been made in the decays H + —>■ cs (and in the charge 
conjugate modes), but this final state is of interest only in the region tan/? < 1, which we 
do not entertain here. However, as already mentioned in the introduction, the characteristic 
polarisation of the r produced in the decays W ± — > t ± v t and H ± — > r ± v T ^ which reflects 
itself in the energy distributions of the r ± -decay products, can be used to discriminate the 
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Fig. 8: Lepton energy spectra from the decays t — > bH + — * b{T + v T ) versus the r energy fraction 
x = x T for M H + = 120 GeV and tan/3 = 22. The solid curve shows the Born spectrum (see (|55|) ) 
and the dashed curve is the spectrum including the (QED + QCD) radiative corrections from 
Fig. [7] (i.e. QCD corrections are taken from Eq. (|B24j) and the QED corrections are taken from 



iy ± -induced and if ^induced final states. In this section, we calculate the energy spectra 
of the so-called single charged-prong events in r-decays. The -^-polarisation effects on the 
r ± decay products have already been investigated in the literature, in particular in 0, 
which we shall make use of, convoluting these spectra with the r + -energy spectrum from 
the decay chain t — > b(W + , H + ) — ► 6(r + z/ r ) calculated by us here. To that end, we consider 
the following r decay chains 

t b(W + , H + )tv t -> jet(b) +u t + p t + Ui + 1, 1 = e + , ; (60) 

t -> b(W + , H + )tv t -> jet(b) + v T + v T + 7T + ; 

t -> b(W + , H + )tv t -> jet(b) + v T + u r + vr + + 7T°; 

t -> b{W + , H + )tv t -> jet(6) + z/ T + z/ r + 7T+ + 2tt ; 

i -> &(W + , # + )rz/ r -> jet(6) + z/ T + z/ r + 2tt + + ?T , 
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involving the leptons e + , the 7r + , the vector and the axial- vector mesons p + and aj", 
respectively, with the subsequent decays of the p + and af, as indicated. Keeping in mind the 
long-distance nature of the QED interactions, providing the " large logarithms" , one must 
include the structure function associated factors only with the final charged particles-leptons 
or pions (see Fig. [2, b and Fig.[TU a, b). In the rest frame of the top quark, the r-leptons from 
the decays t — > b(W + , H + ) — > b{r + v T ) have much larger energy and 3-momentum compared 
to ther- mass, i.e. E T ^> m T . The energy spectrum of the r-lepton decay products must be 
modified to take this into account [36]. For example for the decay r — > pvv of the r-lepton 
with energy E T , the E^-energy spectrum can be obtained from ( fT5l) (see also Eq. (2.8) in 



i 



dBr 1 ^^ f , n f dBr B r ( x 



dcos9 / — — — — S [ z (1 — cos9) ) dx 



dz J J d cos 9dx V 2 

z 

TD 

= $ z ( z ) = —(l-z)[5 + 5z-Az 2 -P T (l + z- 8z 2 )] , (61) 
3 

where z is the energy fraction of the fi in the indicated decay: 

,= ! = £(!-„,, _|, (62) 

x = 2E t j - /m T , and 9 is the angle between the directions of the r and the /i 3-momenta. The 
index a in E a shows the final particle involved. Here, a = fi (the expression above holds 
also for the e + -energy spectrum). We also need the energy spectra for other particles in the 



r-lepton decay, i.e. for a = n + ,p,ai. The corresponding distributions were obtained in [36]: 
for $„. (z) see Eq. (2.4) and for $ p ai (z) see Eq. (2.22) in the cited paper. 



A. Leptons in final state 

For the r + -decay with the leptons in the final state (i.e. a = e + , p + ), the final expression 
for the lepton energy spectrum in the Born approximation is 

dBr^+^ + _ f dy T dBr^^ / xA ^ 



dxt J y T dy T \y 



j-i 



where the t-quark decay width — is taken from (1361) or (|55|) and (z) was defined 

in (|6T|) . To take into account the QED radiative corrections in LLA we again use the SF 
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Fig. 9: Differential branching ratio 



dBr 



t-t b W+ — 6 ( ( t + -. fi + v T ) v T 



dx u 



(see ([63]) ) as a function of the jj^ 



energy fraction in the Born approximation (solid curve) and with the QED and QCD radiative 
corrections (dashed curve) (see (|64l) ). 

approach illustrated in Section [III (see formula (jSJ) and Fig. [21 b). The radiatively corrected 
spectrum is then given by 



dBr, 



QE D LLA +QCD 

dxi 



"pfot 
1 t,RC 



—D, 

y 



Xl 

— )« 
y 



d%h x 



2v^ 



dT 



Born 



a. (2 FW(^ 



dxbdy 

where the first entry in the curly braces is for the decay channels which go via the W + b 
intermediate state, and the second entry is for the decays which go via the intermediate 
state H + b. The function F w (xi,Xb) represents the non-leading contributions of the QCD 
corrections and was defined in fj40|) . The definition of the function Fh(xi,Xi,) is given 
in (1B24j) . The differential branching ratios from the decay chain t — > bW + — > b(r + — > 
^ + v ll v T )v r are shown in Fig. [9] and from the decay chain t — > bH + — > b(r + — > ii r v i Jj T )v T in 
Fig. EH 
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Fig. 10: Differential branching ratio 



dBr 



t -> b H f -> b ( ( t t — > p. "I- V/1 v r ) v T 



as a function of the fi + energy fraction 



in the Born approximation (solid curve) and with the QED and QCD radiative corrections 
(dashed curve), involving the bH + intermediate state. 

B. Hadrons in the final state 

In this case, we have to take into account the decay chain involving the final decays 
t + — > (n + , p + , af)u T with the subsequent decays of the p + - and the a^-mesons into pions. 



The r-energy spectrum is already given above. Consider first the decay r + 
case the p + -energy spectrum is given by 

i 



dBr 



B 



oc 'j- dec j- V 0bq- 



p + z/ T , in which 



(65) 



where the function <& p (x) describes the conversion of the energetic r-lepton into the energetic 
p-meson. This function was calculated in Ref. [36j, which we incorporated in our numerical 
calculations. The distribution in the pion energy fraction x n resulting from the p + — > 7r + 7r° 
decay takes the form: 

(66) 



dBr 



t-*b+u T +(r-*P T +(p-nr+ ... )) 



Born 



dx p dBr Boin ^ / x^ 

oc p doc p \^ X p 
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where the function R p (x) describes the conversion of the energetic p meson into energetic 
pions (i.e. p ± — > 7r ± 7T°): 



R (x) = 



x 



X, 



(67) 



This function was investigated in 37| (see Fig. 1 in 37j]). 

The radiative corrections ("large distance contributions") can be obtained by using the 
structure function approach as: 



dBr 



t-^b+V T + {T^D T + { P ^TT+ ... )) 

RC 



{l,B?} fdy n 



T^tot 

1 t,RC 



X 



dr 



t-^b+V T +(T^i> T +(p-*TT+... )) 

Born 

dy n 



D I — , 



where D n (z) is the structure function of the charged pion 54j : 



D w {x, A,) =S(l-x)+ /3P« (x) + -P 2 P^ ( x ) + 



and the quantities Pi™^ (x) have the form: 



p(i) 



2x 
1 — x 



lim 



2x 

1 — X 



9 (1 - x - A) + (2 In (A) + 2) 5 (1 - x) 



The formula similar to ( 1591) in the case of pions reads as: 

i i 
J ^P« (~)<My) = 0(x){21n(l-x)+2} + | 



2x 
y y-x 



(68) 



(69) 



(70) 



[0(</)-0(x)], (71) 



where 4>(x) is any arbitrary function we want to convolute with the structure function 
D w {x, Pn) in the leading order perturbation theory. The smoothed form of the structure 
function (x, /3 n ) takes the form: 



D n (x,P n ) = 2P w (l-z 



2/3^-1 



(72) 



The energy distribution of the charged pion obtaining from the decay a~l — ► 7r + ... from 
the parent r decay can be derived analogously. 
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6.) 



Fig. 11: Kinematics of the Structure Function method for the decay a) t — > 6z^ T ^ T 7r°7r + and b) 
i — > 6P T i/ T 7r + 7T _ 7r + , involving a 7r + in the final state. 




Fig. 12: Differential branching ratio ^ (see (I63p as a function of the tt + 

energy fraction in the Born approximation (solid curve) and with the QED and QCD radiative 
corrections taken into account (dashed curve) (see (f6lj) ). 

VII. SUMMARY 

In the first part of our paper, we have calculated the QCD and QED radiative corrections 
to the semileptonic decays t — > bW + — > bl + vi [t = e, fi, r) in the SM. Of particular interest 
are the charged lepton energy spectra, which we have calculated using the SF approach to 
resum the leading order (in QED) and leading and next-to-leading order (in QCD) contri- 
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Fig. 13: Differential branching ratio ^ as a function of the ir + energy fraction 

x n in the Born approximation (solid curve), and with the QED and QCD radiative corrections 
taken into account (dashed curve), involving the bH + intermediate state. 

butions. These spectra will be measured accurately at the LHC and will be crucial to check 
the lepton (e, fi, r) universality in the semileptonic decays of the top quarks in SM. In doing 
this, it will be crucial to take into account the QED and QCD radiative corrections in the 
energy spectra. The numerical extent of such corrections is shown in Fig. [5] for the ratios 
R er and R^ T1 which is one of our principal results in this paper. The rest of our paper is 
addressed to the possible effects of a charged Higgs boson H ± with Mjj± < m t — mj in 
the semileptonic decays of the top quark. To avoid the constraints on M H ± coming from 
the B — > X s 7 decay, we assume that the Higgs sector is part of a supersymmetric theory. 
Except for the SUSY radiative corrections, which can be effectively taken into account by 
the supersymmetric renormalisation of the 6-quark mass, there are no other effects of the 
supersymmetric sector on the decay widths and distributions. We have considered only the 
large-tan f3 parameter space of this model, in which case the decays of the H are dominated 
by the final state H ± — > r v T . In Figs. H] and [HI we have contrasted the Born and radiatively 
corrected r-lepton energy spectra from the decays t — > b(W + ,H + ) — > t + u t for a specific 
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Fig. 14: Differential branching ratio as a function of the pion energy 

fraction x w in the Born approximation (solid curves) and including the QED and QCD radiative 
corrections (dashed curves) . The upper curves are for the transverse polarisation of the p and the 
lower curves correspond to the case where the p is longitudinally polarised. 

choice of the parameters M H ± = 120 GeV and tan/3 = 22. While the Born level spectra are 
well documented in the literature, effects of the radiative corrections on the spectra are, to 
the best of our knowledge, new results. 

The contribution of an H ± in t (i) decays, if allowed kinematically, will enhance the decay 
rate for t — > br + u T (t — > 5t~z/ t ), which is the main i? ± -search strategy at the Tevatron. 
However, with a much larger tt cross section and the luminosity anticipated at the LHC, 
this search strategy can be further strengthened by taking into account the different r ± - 
polarisations in the decays W ± — > ^u T and H ± — > r ± z/ r . As the polarisation information 
of the is transmitted to the decay products of the r 1 * 1 , we have calculated the energy 
distributions of the charged particles (e ± , p ± , tt ± , p ± , af) in the single-charge-prong decays of 
the t ± , as well as the inclusive charged pion spectra from the decay chains t — > b(W ± , — > 
6(r ± , v T ) — > b-K ±J i-X. The results at the Born level are well known in the literature. We have 
calculated the perturbative stability of these distributions. The entire effects of the radiative 
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Fig. 15: Differential branching ratio as a function of the pion energy frac- 

tion x w in the Born approximation (solid curves) and including the QED and QCD radiative 
corrections (dashed curves), involving the bH + intermediate state. The upper curves are for the 
longitudinal polarisation of the p and the lower curves correspond to the case where the p is 
transversely polarised. 

corrections presented here can be implemented in existing Monte Carlos, such as PYTHIA 
and HERWIG, to provide an improved theoretical profile of the semileptonic decays of the 
top quark in the SM and can be combined with FeynHiggs to include the SUSY-related 
corrections specific to particular MSSM scenarios. 
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APPENDIX A: NUMERICAL VALUES OF THE INPUT PARAMETERS 



For our numerical calculations we used the following values of the parameters: 
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0.040 
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x 

71 



t-»6W+— i>((aJ--»(37r)+)i/ T ) 

Fig. 16: Differential branching ratio as a function of the pion energy 

fraction in the Born approximation (solid curves) and including the QED and QCD radiative 
corrections (dashed curves). The upper curves are for the transverse polarisation of the a\ and the 
lower curves correspond to the case where the a\ is longitudinally polarised. 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 


a' 1 


137.035999679 


m t 


171.2 GeV 


tan P 


40 




a s 


0.1176 


m h 


4.20 GeV 


Br (r -»■ iiv T v^) 


17.36 x 10" 


-2 


m e 


0.510998910 MeV 


M n 


0.13957018 GeV 


Br (r — > itv T ) 


10.91 x 10" 


-2 


«V 


105.6583668 MeV 


M p 


0.775 GeV 


Br (r -> pv T ) 


25.52 x 10" 


-2 


m T 


1.77684 GeV 


Tp 


0.1462 GeV 


Br (r — > aiU T ) 


1.859 x 10" 


-1 


M w 


80.398 GeV 




1.230 GeV 








Tw 


2.141 GeV 




0.420 GeV 








M H 


120 GeV 


G F 


1.16637 x 10~ 5 GeV" 2 








T H 


2 GeV 


9 


0.653057 
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Fig. 17: Differential branching ratio 



dBr 



:^6if+->i)((aJ"->(37r)+)i/ T ) 



as a function of the pion energy frac- 



tion x w in the Born approximation (solid curves) and including the QED and QCD radiative 
corrections (dashed curves), involving the bH + intermediate state. The polarisation of the a\ is 
indicated in the figure. 

APPENDIX B: RADIATIVE CORRECTIONS TO TOP QUARK DECAY VIA 

CHARGED HIGGS 



Here we give details of the QCD radiative corrections to the width of t-quark decay 
t{p) — > b(pb) + r(p T ) + u{p u ) with the charged Higgs boson in the intermediate state. 

The lowest order QCD corrections can be calculated in a similar way as in QED and in 
the final result one must do the replacements 

TV 2 — 1 4 

a^a s C F , C F = -^- = -, (Bl) 

where N c = 3 is the number of quark colours. 

We start from the counter-terms associated with the t and b quarks. Taking them into 
account yields a multiplicative renormalisation factor in the expression for the differential 
width 



dT — >• dTZ bt} 



(B2) 
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Fig. 18: Normalised inclusive pion energy spectra from the decays t — > b(W + , H + ) — ► 6(7r + + X)P t 
as a function of the pion energy fraction x w in the Born approximation (solid curves) and including 
the QED and QCD radiative corrections (dashed curves). 



m — 



, A 2 3 , m 2 9 , to? 
ln ^f + 2 ln 4 + 2- 21n ^ 



(B3) 



2vr 

where m t and m b are the masses of top and bottom quark. The auxiliary parameters A and A 
are introduced to regularise the ultraviolet (UV) and infrared (IR) singularities, respectively. 
Sometimes A is also dubbed as a fictitious "photon (gluon) mass". The dependence of the 
decay width on these parameters will disappear from the final result. The UV-cutoff A will 
be absorbed by the coupling constant renormalisation and the IR-cutoff A will be cancelled 
by taking into account the emission of the virtual and real gluons. 

Virtual corrections associated with the vertex type Feynman diagram require the calcu- 
lation of the following integral involving the loop 4-momentum 

v= f d 4 k 7^ - k + m h ) (p - k + m t )^ 

J in 2 (k 2 - \ 2 )((p b -ky - m 2)((p-jfc)2 -m 2 )' 1 ' 

Using the Feynman prescription of combining the denominators and performing the loop 
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momentum integration (here we must impose the ultraviolet cut-off) we arrive at 



V = 4 In 



2 2 

— -2(pp b ) [ -f ln^f-2(m t + m b ) / -| (zm 6 + (1 - x)m t ) , (B5) 

fit' Mb J Px A J Px 



with p 2 = x 2 m\ + (1 — x) 2 m 2 + 2{pp b )x{\ — x) and the unit matrix in the Dirac space is 

implied. Below we use the explicit form of the 1-fold integrals 
i 



dx 
P 2 x 



l;(l-x);ln 



Px 



mr 









" \ m ( 




mfy \ 



ym t 



\ny; 



1 



2 1 , 2 m t 
, m y — - m — ^ 

mfy V 2 rrr 



2 ( 1 

2Li 2 1 - - 

6 V y 



mftl-y)-*' ~2.. » ™2 — ' ' (B6) 

The next step consists of the calculation of the contribution arising from emission of real 
gluons - soft and hard ones. Standard calculation for the case of soft gluon emission u < 
AE b <C m t (we work in the rest frame of top quark) leads to 

dT so ft 



dT, 



4tt 2 



J d 3 k 


( p 


_ Pb \ 




[pk~ 


Pbk) 



a 

7T 



2(/-l)ln^ + l + /-/ 2 -^ 
A o 



I = In 



ymt 
m b ' 



(B7) 



Extracting the factor Z = 1 + |^ In mt _ t2 and using the ultraviolet-regularised quantities we 
can write 



ZdT, 



dT. 



(B8) 



Collecting the contributions of the Born level and the virtual and soft real corrections, we 
obtain: 



dT = dT B [l + ^(L-l)(2\nA + h + -[-\nA-^ + l\n 2 y-^\ny + Lt 2 (l--) 

ITT I 71 III y 



7T 



2 1 



■ m y}} 



(B9) 



6 l-y 

withA = AE b /E b ,L = ln^. 

Note that the term containing In A are connected with the emission from the "light" 
6-quark and the heavy top-quark. 

— # 

Let us now consider the emission of the hard gluon with momenta k — (uj,k),k > AE. 
The relevant phase volume 



d$ 4 



(2tt) 4 d 3 p b d 3 p T d 3 p v d 3 k 
(2tt) 12 2E b 2E T 2E U 2uj 



S A (P-Pb-Pv -Pr~ k), 
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can be transformed as 
m 



4 



2 13 

where 



-—^-dyx T dx T dxdzdO T 5 ( I — x T — x — y + z -\ — — j~ ^ ) , (BIO) 



and 



,^ 2dcidc 2 ,„ \ 

dO T = -., (B12) 

a/1 — c 2 — c 2 — c 2 + 2cic 2 c 

is the angular phase volume of the r lepton, and 

c = cos(k,p b ), ci = cos(k,p T ), c 2 = cos(p T , p b ). (B13) 
Explicit calculation yields 

[ dO T 5 (l- Xr - x - y + z+ 2p ^ k + Pb) ) = iL , = ^ + y )2 _ 4z. (B14) 

y V m i J x tR 

and the variable z is bounded by 



z m < z < xy, z m = 1 < 1. (B15) 

ymf 

Summed over the final spin states, the matrix element squared leads to 

xy 

dT uncoll a f dx f dz dT(x + y-z) 1 ^ + ^ + ^ ^ 



dx T dy 2n J J dx T dy (x + y — z) R 

where 



F 2 



I/ 2 + + a; ) 2 2m 2 x + y 



2 JL 

x A 



xz m 2 ^ 2 



= — (l + y + x) + 4( 2 + ^) ■ 

a; ar 

It is convenient to introduce the small auxiliary parameter a (z m <C cr <C xy ~ 1) and 

— * 

extract the contribution of the collinear kinematics k\\pb- 

Only Fi gives the contribution in the collinear region (z < a): 

(BIT) 



a f dt 
dT coll = — I jdT (t) 

y(l+A) 



1 - f — - ■ - t 
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with 



In 



m t ya 



xra 



2 ' 



A 



2AE 

m t y ' 



Contribution of F± from the non-collinear region can be put in the form: 



a 
2^ 



dt 



dT (t) 



y(l+A) 



y{t-y) 



I"! 



a 



dz 



£(*,*)] 



(B18) 



with 



E(t,z) 



dF(t - z)t 2 



1. 



(t - z)dY(t)^t 2 - Az 

Note that the second term containing £(£, z) is finite in the limit mj 
The contribution of the second term (F 2 ) can be cast in the form: 

yit-y) 



(B19) 



0. 



a 

7T 



dt y 2 dT(t) n 
/ // £ 2 rr 



J/(l+A) 3/ U 

The first term above combined with the term — -dr(y) In A (see (IB9I) ) gives a quantity which 
is finite in the limit A — > 0. Emission from the light quark has a form predicted by the 
Structure Function approach. Combining all the contributions, we arrive at the following 
expression for the QCD corrected double (Dalitz) distribution in the variables x T , y: 

i 

>if /» - n dT(t,x T ) f i , a s ^ } 



dT (y,x T ) 



dx T dy 



dtdx T 



l + —F H 

7T 



where Fjj is the .KT-factor which contains all the non-enhanced terms. On integrating the 

6-quark energy fraction, the mass singularities ($(y) = f^(lia(y 2 m 2 /mf) — 1)) for m b — ► 

will disappears due to the relation 

11 i 
JdyjjD (|,/?) F{t) = J F(t)dt. (B22) 

y 

Thus, in an experimental setup involving an averaging over the 6-jet production, the resulting 
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r-meson energy spectrum fraction is described by the following expression: 



dr 

dx T 
dT 



dxh 



dT 



dxbdx T 



l—x T 



dT 



B 



OL s Gt, 



dtdx T dtdx T 



71 



, dT B y 2 1 
dy —, — h 



dydx T t 2 t — y 



+ 



ol s Cf dT b 

71 dtdx T 

dT b 



*(1+A) 

h - In 2 t - - In t 

2 2 2 



lnt 
1 - t 



C 2 + Li 2 ( 1- - ) - --InA 



t 



dtdx 



- (l - —F H (x T ,t) 

T V 71 



(B23) 



As anticipated, this expression does not depend on the small auxiliary parameter A C 1. 
Thus function F H (x T ,Xb) is defined as 



hit 



F H (x T ,t) = C F { - - - ln 2 t+ -lnt + + ( 2 - Li 2 1 - - + - + lnA 



+ 



(■ 



dT 



2 2 

-l 

B 



V dtdx T 



X 



t(l+A) 



1 -t 

dr B y 2 i 

dydx T t 2 t — y 



1 



t 



(B24) 
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